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A  DIPOLE  APPROXIMATION  OF  THE  BACKSCATTERING 
FROM  A  CONDUCTOR  IN  A  SEMI-INFINITE  DISSIPATIVE  MEDIUM 
WITH  APPLICATION  TO  SUBMARINE  DETECTION 


Prepared  by: 
Martin  B.  Kraichman 


ABSTRACT:  The  backseat  ter  Ing  of  a  uniform  plane  wave  by  a 
conductor  In  a  seml-lnfinlte  dissipative  meditim  is  discussed. 

The  conductor  is  assumed  to  act  as  both  an  electric  and  a 
magnetic  dipole  with  moments  which  are  obtained  from  the 
electric  magnetic  polarisabllities  of  the  conductor,  respectively. 
Using  these  induced  moments,  expressions  are  derived  for  the 
backscattered  electric  field  at  a  point  on  the  surface  of  the 
dissipative  half-space  directly  above  the  dipoles.  Both  harmonic 
and  transient  excitation  are  considered.  Numerical  results 
are  presented  for  the  backseat te'^ed  electric  field  from  a  per¬ 
fectly  conducting  finite  cylinder  in  sea  water  when  the  cylinder 
is  excited  by  a  single  sferic  pulse.  The  application  of  these 
results  to  the  detection  of  suomarines  is  discussed. 
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INTRODUCTION 

1.  The  scattering  of  electromagnetic  waves  by  conductors 
In  an  Infinite  homogeneous  dissipative  medium  has  been  treated 
by  various  investigators  in  recent  years.  The  classical  prob¬ 
lem  of  the  scattering  of  a  plane  wave  by  a  spherical  object  is 
reviewed  by  Stratton  (1941).  Wait  (1951,  1960)  has  considered 
the  response  of  a  con^ctlng  sphere  to  a  uniform  time  varying 
magnetic  field  and  to  the  fields  of  electric  and  magnetic 
dipoles.  The  scattering  by  an  infinite  inhomogeneous  conducting 
cylinder  under  the  influence  of  s  time  varying  magnetic  field 
has  been  investigated  by  Negi  (1962). 

2.  In  many  practical  situations,  it  is  of  interest  to 
consider  the  scattering  of  electromagnetic  waves  by  conductors 
in  a  semi-infinite  homogeneous  dissipative  medium.  Recently, 
Galejs  (1962),  using  a  dipole  approximation,  treated  the  prob¬ 
lem  of  the  scattering  from  a  conducting  sphere  in  such  a  medium 
when  the  sphere  is  excited  by  a  surface  wave  or  by  fields  from 
vertical  electric  or  horisontal  magnetic  dipoles  in  the  lossless 
half -space.  The  analysis,  however,  is  restricted  to  field 
points  with  horizontal  ranges  from  the  sphere  which  are  much 
greater  than  a  wavelength  in  the  dissipative  medium. 

3.  The  present  paper  considers  the  scattering  of  a 
uniform  plane  wave  by  a  conductor  of  finite  dimensions  embedded 
in  a  semi-infinite  dissipative  medium.  The  frequency  of  exci¬ 
tation  and  the  conductivity  of  the  medium  are  such  that  the 
displacement  current  can  be  neglected.  A  dipole  approximation, 
in  which  the  electric  and  magnetic  dipole  sMxnents  of  the  con¬ 
ductor  are  obtained  from  the  electric  and  magnetic  polariza¬ 
bilities,  respectively,  is  used  to  avaluate  the  backscattered 
electric  field  in  the  interface  directly  above  the  conductor. 

The  electric  field  step  response  is  obtained  from  the  response 
to  harmonic  excitation  by  transforming  the  Fourier  integral 
into  a  convolution  integral.  Graphs  of  the  electric  field 
step  responses  are  presented  for  the  backseat ter ing  from  both 
electric  and  mametic  dipoles.  Numerical  results  are  pre¬ 
sented  for  the  backscattered  electric  field  from  a  perfectly 
conAieting  finite  cylinder  in  sea  water  when  the  cylinder  is 
excited  by  a  single  sferic  pulse.  The  application  of  these 
results  to  the  detection  of  submarines  is  discussed. 
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DIPOLE  APPROXIMATION 

4.  A  uniform  plane  wave  le  aaeumed  to  propagate  vertically 
In  a  semi-infinite  dissipative  medium  having  a  norlzontal  planar 
Interface.  The  displacement  current  In  the  dissipative  medium 
Is  considered  to  be  negligible  compared  with  the  conduction 
current.  The  upper  half-space  Is  assumed  to  have  the  properties 
of  free  space.  A  perfectly  conducting  body  (or  one  with  negli¬ 
gible  field  penetration) ,  which  Is  embedded  In  the  dissipative 
half -space,  will  have  a  surface  current  Induced  on  It  by  the 
Incident  plane  wave.  The  resulting  scattered  electromagnetic 
field  may  be  approximated  by  that  crom  a  horizontal  electric 
dipole  and  a  horizontal  magnetic  dipole  if  the  vertical 
dimension  of  the  conductor  Is  much  smaller  than  a  wavelength  In 
the  dissipative  medium  and  If  the  dimensions  of  the  conductor 
are  small  compared  with  the  distance  to  a  field  point.  In 
this  approximation,  the  electric  field  of  the  plane  wave  induces 
a  horizontal  electric  dipole  moment  In  the  conductor  while  the 
magnetic  field  Induces  a  horizontal  magnetic  dipole  moment. 

Both  of  these  moments  act  as  a  source  of  electromagnetic  field 
which  can  be  determined  from  the  propagation  equations  for 
horizontal  electric  and  magnetic  dipoles  in  a  semi-lnflnlte 
conducting  medium. 


Induced  Electric  Dipole  Moment 


5.  Referring  to  the  rectangular  coordinates  shown  In 
Figure  1,  the  surface  of  a  semi -Infinite  dissipative  medium 
Is  located  at  z  -  o.  A  uniform  plane  wave,  with  harmonic 
variation  eJ"^,  propagates  In  the  medium  in  the  positive  z 
direction.  The  electric  and  magnetic  fields  are  given  by 
Ex(*><^)  *od  Hy(z,w)  respectively.  If  the  Incident  electric 
field  Is  substantially  uniform  over  the  height  (z-diaenslon) 
of  a  perfect  conductor  at  a  depth  z  «  d,  the  conductor  will 
acquire  an  Induced  electric  dipole  moment,  Px(d,ti)),  which  may 
be  expressed  as 


^0 


where  (  ■  permittivity  of  the  dissipative  medium 

9  m  conductivity  of  the  dissipative  medium 

■  electric  polarizability  of  the  perfectly 
conducting  body 
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Incident  electric  field  at  a  depth  d. 

The  quantity  enclosed  In  brackets  represents  the  complex 
permittivity  of  the  dissipative  medium.  If  the  displacement 
current  In  the  dissipative  medium  Is  neglected  (  <»■  ), 

then  (1)  reduces  to 


6.  The  current  moment  P»(d,<A))  may  be  obtained  by  multi¬ 
plying  the  right  side  of  (2)  by  J  w  .  Thus, 

=  (ra*  Ex  (3) 


and  the  Induced  electric  current  moment  Is  seen  to  be  In  phase 
with  the  Incident  electric  field.  Multiple  reflections  between 
the  conducting  body  and  the  surface  of  the  half-space  will  be 
considered  later  In  this  section. 


Induced  Magnetic  Dipole  Moment 

7.  The  permeability  of  the  dissipative  medium  and  the 
conducting  body  is  assumed  to  be  that  of  free  space.  If  the 
incident  magnetic  field  is  substantially  uniform  over  the 
height  of  the  conductor  at  a  depth  d,  then  the  conductor  will 
acquire  an  Induced  magnetic  dipole  moment,  My(d,w),  which  Is 
opposite  In  direction  to  the  Inducing  magnetic  field.  The 
dipole  moment  Is  expressed  as 


where  magnetic  polarizability 

Hy(d,w)-  Incident  magnetic  field  at  depth  d. 

Since  My(d,(«))  is  opposite  In  direction  to  Hy(d,w),  the  polariza¬ 
bility  will  be  negative.  The  magnetic  current  moment, 
l4(d,«J),  Is  defined  as 
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cn 


where  -  permeability  of  free  space. 

8.  The  polarizabilities,  and  ,  are  functions  of 
the  volume  and  shape  of  the  conductor  and  the  relative  orien¬ 
tation  of  the  conductor  with  the  inducing  field.  The  expres¬ 
sions  for  the  current  moments  given  by  (3)  and  (5)  agree  with 
those  given  by  Galejs  (1962)  for  a  conducting  sphere  of  volume 
V  if  the  well  known  values  or^  >  3v  and  -3/2  v  are  used  for 

the  polarizabilities. 


Perturbing  Effect  of  Reflections 

9.  The  incident  plane  wave  at  the  conductor  will  be 
perturbed  by  reflections  of  the  backscattered  electric  and 
magnetic  fields  from  the  surface  of  the  dissipative  half-space 
The  perturbing  effect  of  a  single  reflection  on  the  electric 
and  magnetic  current  moments  oz  spheres  and  certain  cylinders 
will  be  approximated  and  shown  to  be  small. 

10.  After  a  single  reflection  from  the  surface  of  the 
dissipative  medium,  the  reflected  electric  and  magnetic  fields 
Exr(<l>^)  and  Hyr(d,t<)),  respectively,  induce  electric  and 
magnetic  current  moments  in  the  conductor.  The  moments  are 
given  by 


(i) 

and 

Mjr  =  j**y**®t*n  Hijr  (d,uf)  ,  (?) 
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The  reflected  electric  end  magnetic  fields  In  (6)  and  (7)  may 
^  expressed  In  terms  of  the  current  moments  F^xl(<^i*^) 

^1  (<!>><))  which  are  Induced  In  the  conductor  by  the  Incident 
plane  wave.  Because  of  the  large  Impedance  mismatch  at  the 
surface  of  the  dissipative  medium,  It  Is  assumed  tluit  the  back- 
scattered  electric  field  Is  totally  reflected  In  opposite 
phase.  In  returning  to  the  conductor  at  a  depth  d,  the 
reflected  wave  will  have  experienced  a  spherical  spreading 

distance  2d.  Using  these  approximations,  the  reflected 
fields  may  be  obtained  from  the  well  known  propagation  equations 
for  the  electromagnetic  fields  from  electric  and  magnetic 
dipoles  In  an  Infinite  dissipative  medium.  The  electric  dipole 
contributes  to  both  the  reflected  electric  and  magnetic  fields. 
Similarly,  the  magnetic  dipole  also  contributes  to  both  fields. 
The  reflected  fields  are  therefore  written 


<f) 


and 


Hjjr  r 


R«l 


e 


where  r  -  (J <»>/*,  o*)^. 

For  |Tdi»  1,  the  exponential  attenuation  factors  in  (8)  and  (9) 
will  ensure  the  smallness  of  the  perturbing  electric  and  mag- 
netlc  fields.  For|r^|«  1,  the  perturbing  moments  in  (6)  and 
(7)  may  be  expressed  as 
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and 

Mjr  (d.U>) 
(lu)) 


'  (««+  » 


where  use  la  made  of  the  expressions  In  (8)  and  (9)  and  of 
the  definitions  in  (3)  and  (5). 


11.  The  ratios  in  (10)  and  (11)  may  readily  by  shown  to 
be  small  for  spherical  conductors  and  for  cylinders  which  are 
neither  extremely  needle  hor  extremely  disc  shaped.  For  a 
sphere  of  volume  v  and  radius  a,  oc^  ■  3v  ■  4  r  a3  and  • 

The  ratios  in  (10)  and  (11)  then  become 


PkI  (<I,w) 

and 

H'ji- 

M',-.  (iM) 


-  f  ^  (i-rd) 

■  t  T’  • 


Since  it  has  been  assumed  that  a«  d  and  that  |r4|<ICl>  it 
follows  that  the  perturbations  in  the  electric  and  magnetic 
moments  are  small  for  a  single  reflection. 

12.  The  electric  and  magnetic  polarisabilities  of  a 
right  circular  cylinder  are  gi.ven  in  Table  1.  The  subscript  / 
denotes  longitudinal  excitation,  i.e.,  with  the  particular 
field  vector  of  Interest  along  the  axis  of  the  cylinder.  The 
subscript  t  denotes  an  excitation  which  is  transverse  to  the 
axis  of  the  cylinder.  The  values  in  Table  1  are  obtained  from 
Taylor  (1960a,  1960b).  The  entries  for  a/b  -  0.1  represent 
values  which  ^ve  been  extrapolated  from  the  results  in  the 
above  references.  In  the  dipole  approximation,  the  dimensions 
of  the  cylinder  must  be  much  smaller  than  the  depth  at  which 
the  cylinder  is  embedded,  so  that  v«d^.  It  is  therefore 
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evident  from  (10),  (11),  and  Table  1,  that  the  moments  due  to 
the  singly  reflected  electromagnetic  field  are  small  compared 
with  the  moments  due  to  the  Incident  plane  wave  for  cylinders 
which  are  neither  extremely  needle  nor  extremely  disc  shaped. 

13.  Higher  order  perturbations,  resulting  from  multiple 
reflections,  will  be  even  smaller  than  the  perturbation  due  to 
a  single  rejection. 


BACKSCATTERED  ELECTRIC  FIELD  AT  THE  SURFACE 

14.  In  this  section,  expressions  will  be  derived  for  the 
backscattered  electric  field  at  the  surface  of  the  dissipative 
half -space  for  a  field  point  directly  above  the  conductor. 

The  electric  field  contributions  from  the  Induced  electric 
and  mametlc  dipoles  will  be  treated  separately  and  both 
harmonic  and  transient  excitation  will  be  considered. 


Electric  Field  From  an  Induced  Electric  Dipole 
Harmonic  Excitation 


IS.  Von  Aulock  (1952)  has  derived  expressions  for  the 
electric  field  In  a  seml-lnflnlte  dissipative  medium  In  a  region 
directly  above  an  embedded  horlsontal  electric  dipole.  As  Is 
shown  In  Figure  1,  the  dipole  Is  situated  at  the  point  (o,o,d) 
and  Is  directed  along  the  X-axls.  The  dipole  Is  excited  at  an 
angular  frequency  w .  At  the  coordinate  origin  In  the  surface, 
the  electric  field  E**  (o,w),  Is  directed  along  the  X-axls  and 

Is  written 


iirrtl' 


where  ir>  (Jw/i,9’)^  and  K2(l^d)  Is  the  modified  Bessel  function 
of  the  second  kind.  Using  (3)  and  the  plane  «Mive  relationship 

H,(o.4))  (liT) 

which  relates  the  electric  field  Incident  at  the  conductor  to 
the  magnetic  field  at  the  surface  of  the  half-space,  (14)  may 
be  written 
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Transient  Excitation 


.  16.  The  backscattered  electric  field  In  the  time  domain  , 

Ex  (o>t),  may  be  expressed  as  a  Fourier  Integral  and  written 

.  O') 

The  Fourier  Integral  may  then  be  transformed  Into  a  convolution 
Integral  Involving  the  backscattered  electric  field  response 
at  the  surface  of  the  half-space,  directly  above  the  dipole, 
when  a  unit  step  of  magnetic  field  Is  applied  at  the  surface 
In  the  y  direction.  The  steps  of  the  transformation  are  given 
In  Appendix  A.  Equation  (17)  then  reads 


where  A^  (o,t)  «  backscattered  electric  field  step  response  at 

the  surface  from  an  Induced  electric  dipole, 

Hy  (o,t)  •  time  derivative  of  the  transient  magnetic  field 
^  ^  applied  at  the  surface . 

*(See  Insertion,  page  8a) 

The  step  response  may  explicitly  be  written 


where 


t 


and 


A4(v)  +  SinV 


€  dV. 
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*Slnce  (17)  Includes  all  frequencies,  the  condition  that  the 
displacement  current  be  negligible  compared  with  the  conduction 
current  will  no  longer  hold  for  the  higher  frequencies. 

However,  In  a  medium  with  appreciable  conductivity,  the 
electromagnetic  fields  at  these  high  frequencies  will  be  so 
severely  attenuated  that  their  effect  will  be  negligible  at 
all  but  the  most  shallow  depths.  Since  the  approxl^tlon  In 
(17)  to  the  true  E|*(o,u>)  Is  also  negligible  at  high  frequencies, 
the  results  obtained  by  using  (18)  will  be  of  practical  value. 
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The  above  expression  for  1(|S)  Is  convenient,  for  large  values 
of /I  .  With  further  manipulation,  It  can  be  Ahown  that  I(/9)  may 
also  be  written 


This  form  Is  convenient  for  small  values  of/f  .  The  leading 
terms  of  an  asymptotic  expansion  for  small  a  obtainable  from 
(21)  are  ' 

(1  +  1.68750/9  +  0.95508/*^  +  0.16296/J^ 

-  2.24885/3^  +  .  .  .).  (li) 


,nd  22) .  It  la  evident 


From  the  exponential  factor  e'J  In  (19)-|nd  22) , 
that  the  backscattered  electric  field,  Ex  (o»t).  Is  quite 
small  for  t«  djMo<r.  For  large /S  ,  It  can  be  shown  that  A$  (o,t) 
falls  off  as  t'^..  A  graph  of  the  dimensionless  quantity 
-aJ  (o.t).  4+  9  d^  Is  shown  in  Figure  2.  The  Integrals  were 

evaluated  on  the  IBH  7090  using  Gaussian  Integration. 


Electric  Field  From  an  Induced  Magnetic  Dipole 
Harmonic  Excitation 


17.  Von  Aulock  (1952)  has  also  derived  expressions  for 
the  electric  field  directly  above  a  horlsontal  magnetic  dipole 
In  a  seml-lnflnlte  dissipative  medium.  As  Is  shown  In  Fl^re  1, 
the  magnetic  dipole  Is  situated  at  the  point  (o,o,d)  and  Is 
directed  along  the  Y-axls.  The  dipole  Is  excited  at  angular 
frequency  w .  At  the  coordinate  origin  In  the  surface,  the 
electric  field,  E^*"*  (o,w)  Is  directed  along  the  X-axls  and  is 
written 
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where  /-  And  K]^(/d),  K^Cl^d)  are  the  modified 

Bessel  functions  of  the  second  kind  of  the  first  second  order, 
respectively.  Using  (5)  and  the  plane  wave  relationship 


-r^ 


equation  (23)  may  be  written 

£>,«) •  *  3  <Urd) 


Transient  Excitation 


.  18.  The  backscattered  electric  field  in  the  time  domain, 

Ev  (o,t),  may  be  expressed  as  a  Fourier  Integral  similar  to 
that  In  (17).  The  Fourier  Integral  may  then  be  transformed 
Into  a  convolution  Integral  of  the  form  given  In  (18) ,  where 
AS  (o,t)  Is  now  the  backscattered  electric  field  step  response 
at  the  surface  from  an  Induced  magnetic  dipole.  The  steps  of 
the  transformation  are  given  In  Appendix  B.  The  step  response 
may  be  written  explicitly  as 


\  ^  (I 

erf 


where 


erf(yc)  a  J  ^3 


and 


00 

fly' 

dlfy 

(>•7) 

€ 

(’■*) 
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The  above  expressions  for  Ii(a)  and  IzCfi)  az-e  convenient  for 
large  values  of  A  .  With  further  manipulation,  It  can  be  shown 
that 


where  erfc(x)  ■  /  e  ^  dy.  These  forms  are  convenient  for 

small  values  of/iv  The  leading  terms  of  the  asymptotic  expansions 
obtainable  from  (29)  and  (30)  are 


X/|9)r-l+ (  1  +  1.93750/3  +  0.43945/3^ 

' '  2^(i  3  4  . 

‘  -  0.028687(3'’  +  1.97677(3  +  .  .  .)  ^i) 


and 

e'*  2 

(  1  +  0.18750|3  -  0.076172(>^ 

^  +  0.10071(3^  -  0.23011(3^  +  .  .  .). 

It  can  be  shown  that  (o,tl  behaves  as  e^jv  for  small /3  and 
Is  thus  very  small  for  t«  d^^.o- .  For. large  /3  ,  the  step 
response  can  be  shown  to  fall  off  as  A  graph  of  the 

dimensionless  quantity  -A#  (o,t)  ,  4Tr  o*  d*»  Is  shown  In 

o»  m 

Figure  3.  The  Integrals  were  evaluated  on  the  IBM  7090  computer 
using  Gaussian  Integration. 
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RELATIVE  IMPORTANCE  OF  THE  ELECTRIC  FIELD 
Contributions  From  Electric  and  Magnetic  Dipoles 

19.  The  relative  Importance  of  the  contributions  by  the 
Induced  electric  and  magnetic  dipoles  to  the  backscattered 
electric  field  at  the  surface  will  depend  on  the  relative  magnl 
tudes  of  the  electric  and  magnetic  polarizabilities  of  the 
conductor.  An  estimate  of  the  relative  Importance  may  be 
obtained  by  determining  the  ratio  of  the  Incident  electric 
field  at  the  surface  from  the  magnetic  dipole  to  that  from  the 
electric  dipole.  The  Incident  electric  field  from  the  magnetic 
dipole,  E^  (o,u)),  Is 


and  that  from  the  electric  dipole,  E^  (o,*«)).  Is 

.  (sh) 


Using  the  definitions  In  (6)  and  (7),  the  ratio  of  the  electric 
field  contributions  Is  written 


1  1 't'rci-hY^d' 

When  1,  the  ratio  Is  small,  provided 


From  (3^) ,  It  Is  evident  that 


e;  (v) 

V 

0(g 
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The  contribution  of  the  magnetic  dlpo^  to  the  electric  field 
will,  therefore,  also  be  small  when  j.^  |  !•  This  Is 

the  case  for  a  cylinder  with  a  small  diameter  to  length  ratio 
when  excited  longitudinally  by  an  electric  field.  On  the  other 
hand,  when  | j  Is  sufficiently  large,  the  contribution  of 

the  magnetic  dipole  will  be  larger  than  that  of  the  electric 
dipole.  This  Is  the  case,  for  example,  for  a  cylinder  with  a 
large  diameter  to  length  ratio  when  excited  longitudinally  by 
a  magnetic  field. 


THE  ELECTRIC  FIELD  SCATTERED  BY  A  SUBMARIIC 

20.  It  Is  of  Interest  to  apply  the  results  of  the  previous 
sections  to  the  scattering  of  an  electromagnetic  field  by  a 
submarine  In  sea  water,  since  an  electromagnetic  field  In  sea 
water  Is  exponentially  attenuated  with  the  frequency,  the 
exciting  field  must  contain  sufficiently  low  frequencies  In 
order  to  appreciably  Illuminate  a  deeply  submerged  conductor. 
Man-made  sources  vdilch  radiate  efficiently  In  the  ELF  range  are 
not  at  present  very  practical.  However,  the  ELF  energy  radiated 
by  ll^tnlng  discharges  and  propagated  as  sferics  could  con¬ 
ceivably  be  used. 


21.  A  typical  sferic,  consisting  of  a  quasi-sinusoldal 
portion  followed  by  a  slow  tail,  Is  applied  to  the  sea -air 
Interface  and  Is  assumed  to  pro^gate  vertically  into  the  sea 
water  as  a  uniform  plane  wave.  This  wave  Is  then  backscattered 
by  a  submarine  at  a  depth  of  150  meters.  The  submarine  Is 
approximated  by  a  conducting  cylinder  with  a  diameter  to  length 
ratio  of  1/10  and  a  volume  of  10^  m3.  The  electric  field  of 


the  Incident  tMve  Is  taken  to  be  along  the  longltudl 
of  the  cylinder.  Since  for  this  type  of  excitation 


only,  the  electric  field  contribution  from  the  Induced  electric 
dipole  is  considered.  Although  neglected  here,  the  Induced 
ferromagnetic  moment  of  a  submarine  will,  In  fact,  oppose  the 
Induced  diamagnetic  moment,  thus  helping  to  make  the  contribution 
of  the  latter  moment  even  smaller.  While  it  is  true  that  the 
dimensions  of  the  submarine  are  not  much  smaller  than  the  depth 
of  submergence,  the  dipole  approximation  may  be  expected  to  give 
an  order  of  ma^ltude  estimate  of  the  backscattered  electric 
field  in  the  Interface  directly  above  the  submarine. 


22.  No  attempt  Is  made  here  to  suggest  a  method  for 
detecting  the  backscattered  electric  field.  An  example  of  such 
a  method  Is  discussed  by  Modavls  (1962) . 
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Specification  of  the  Sferlc  Waveform  and  4 

the  Parameters  of  the  Scattering  Model 

23.  The  choice  of  an  analytical  representation  of  the 
waveform  of  the  applied  sferlc  to  be  used  in  the  numerical 
computations  Is  based  on  the  desire  to  approximate  commonly 
observed  waveforms  by  the  simplest  expression  possible. 

Accordingly,  the  oscillatory  portion  Is  given  by 


Hy(o,t)  -  0  for  t  <0  , 

-  Hoe“*^sin(wot-bt2)  for  o  <  t  <  t^, 
.  0  for  t  >  t^ 


where  Hq  ■  3  x  10"^  amp/m 
a  -  2.9  X  10^  sec“^ 
b  »  1.8  X  10^  sec"^ 

a  5.2  X  10^  sec"^ 

-4 

tj.  -  8*6  X  10  sec. 

A  plot  of  the  waveform  specified  In  (38)  is  shown  in  Figure  4 
for  times  greater  than  10"*^  sec.  The  slow  tall  following  the 
oscillatory  portion  of  the  sferlc  Is  represented  as  a  positive 
half-cycle  consisting  of  two  straight  line  portions  Joined  by 
a  section  of  a  parabola.  The  duration  of  the  slow  tall  Is 
three  times  that  of  the  oscillatory  portion  and  the  peak  ampli¬ 
tude  Is  0.4  that  of  Che  maximum  oscillatory  amplitude.  Thus, 


Hy(o,t)  a  «(t  -  tc)  for  tg  <  t  <  t^, 

a  h  - d (t  -  t2)^  for  tj  <  t  <  t3, 

a  -  «  (t  -  tf)  for  t3  <  t  <  tf, 

a  0  for  t  >  t£, 

where  t^  «  2.3  t^ 

^2  “  2.5  tg 

t3  a  2  t2  -  ti 
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-  2  t2  -  tc 

-  ■■  ■  ■ 

0( 

-  0.4  Ho 

A  plot  of  the  above  slow  tall  Is  shown  in  Figure  5. 

24.  The  values  of  the  other  parameters  are  chosen  as: 

<r  ,  conductivity  of  sea  4  mho/m 

Mo  ,  permeability  of  sea  water  >  411  x  10~^  henrles/m 

d  ,  depth  of  scattering  model  150  meters 

P  ,  diameter  to  length  ratio  of  conducting  cylinder 
^  -  1/10 

V  ,  volume  of  conducting  cylinder  •  10^  m3 
Oi  e  I  *  ^ong|tu^lnal  electric  polarizability  >  60 


tf 

(3 

h 


Numerical  Results 

25.  Using  (18)  and  the  scattering  model  described  in 
the  previous  section,  numerical  results  have  been  computed  on 
the  IBM  7090  for  the  backscattered  electric  field  at  the  sea 
surface.  The  background,  against  which  this  backscattered 
electric  field  must  be  detected,  consists  of  the  Incident 
electric  field  at  the  sea  surface.  This  field  may  be  determined 
by  an  application  of  (18) .  The  Incident  electric  field  step 
response  may  be  obtained  from  the  well  known  transmission  line 
equations  for  an  Infinite,  leakage  free,  non-inductive  cable 
in  which  the  voltage  corresponds  to  the  magnetic  field  and 
the  current  to  the  electric  field.  For  example,  see  Goldman 
(1949).  The  result  for  the  Incident  electric  field  step 
response  at  the  sea  surface,  (o,t),  Is  written 
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26.  Graphs  of  the  backscattered  electric  field  step 
response  A®  (o,t),  the  incident  electee  field,  E,^  ^  (o,t), 
and  the  backscattered  electric  field  Ex^  (o,t),  arA  presented 
in  Figures  6,  7,  and  8,  respectively.  Because  of  the  expo¬ 
nential  fall  of  the  step  response,  A|  (o,t),  for  small  values 
of  t,  the  backscattered  electric  Held  at  the  sea  surface  for 
these  times  is  negligible  compared  with  the  incident  electric 
field  at  the  surface.  For  this  reason,  numerical  results  for 
t  <  10**^  sec.  are  not  presented  for  the  backscattered  electric 
field  in  Figure  8.  The  behavior  of  the  incident  Held  and  the 
backscattered  field  for  large  t  may  be  determined  analytically 
from  the  convolution  integral  foraa  for  these  quantities. 

Both  are  found  to  fall  oft  as  for  large  t.  In  Figure  8, 

the  point  marked  ~  3%  of  E,^  4  (0|t)  indicates  the  time  at 
which  the  backscattered  field  is  the  maximum  percentage 
(approximately  3%)  of  the  incident  field. 


CONCLUSION 

27.  If  submarines  at  a  depth  of  500  feet  are  to  be 
detected  in  the  example  Just  given,  a  backscattered  electric 
field  of  the  order  of  lO'^  volts/meter  must  be  measured  at 
the  sea  surface  directly  above  the  submarine.  Obviously,  if 
the  detector  at  the  surface  is  much  off  to  the  side,  the 
value  of  the  scattered  field  will  be  even  smaller.  The  back- 
scattered  field  at  the  sea  surface  directly  above  the  sub¬ 
marine  amounts  at  most  to  about  3%  of  the  incident  electric 
field  at  the  surface.  Such  requirements  are  rather  severe 
and  the  detection  of  submarines  in  the  case  of  single  sferlc 
excitation  therefore  appears  to  be  quite  difficult.  For 
submarines  at  shallower  depths  and  for  different  sferlc 
excitations,  the  detection  problem  should  be  investigated 
further. 
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Table  1 

Electric  end  Magnetic  Polarizabilities  of  a  Right  Circular 
Cylinder  of  Diameter  a,  Length  b,  and  Volume  v  -IT  a^b 

V 


a 

0<«.i 

F 

V 

✓ 

V 

0 

oe 

2.00 

-1.00 

-2.00 

1 

60.00 

2.13 

-1.06 

-1.94 

TIT 

1 

15.1 

2.32 

-1.16 

-1.85 

1 

7 

7.10 

2.61 

-’..31 

-1.74 

1 

3.86 

3.17 

-1.59 

-1.58 

2 

2.43 

4.22 

-2.11 

-1.41 

4 

1.75 

6.18 

-3.09 

-1.27 

oo 

1.00 

-  oo 

-1.00 
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FIG. I  A  CONDUCTOR  EMBEDDED  IN  A 

SEMI-INFINITE  DISSIPATIVE  MEDIUM 
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F1G.2  THE  BACKSCATTERED  ELECTRIC  FIELD  STEP  RESPONSE 
AT  the  surface  from  AN  INDUCED  ELECTRIC  DIPOLE 
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FIG.3  THE  BACKSCATTERED  ELECTRIC  FIELD  STEP  RESPONSE  AT  THE 
SURFACE  FROM  AN  INDUCED  MAGNETIC  DIPOLE 
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FIG.  4  WAVEFORM  OF  THE  OSCILLATORY  PORTION  OF 
THE  APPLIED  SFERIC 
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FIG. 5  WAVEFORM  OF  THE  APPLIED  SLOW  TAIL  SFERIC 
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FIG.  7  INCIDENT  ELECTRIC  FIELD 
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FIG.  8  BACKSCATTERED  ELECTRIC  FIELD 
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FI6.9  CONTOUR  IN  THE  Z-PLANE 
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APPENDIX  A 

1.  By  substituting 

*  jt;  e"““^dx  (m) 

0  ^ 

In  (16),  ths  Inverse  Fourier  trsnsform  In  (17)  mey,  upon 
changing  the  order  of  Integration,  be  written 


=  ["a*  H'  (t^x)  dA  , 


where 


t  •*  (J<t^/U«^)^t  and  K9(Fd)  Is  the  modified  Bessel  function 
of  the  second  kind.  For  convergence  erg  t  In;  >  ±  JT  . 
Thus,  2 


Ml 

e  dij. 


H 


When  t>o,  Win  (A4)  Is  replaced  by  the  complex  variable 
z  >U  Ju  ,  and  the  Integration  Is  performed  around  the  contour 
shown  In  Figure  9 .  When  f o  and  R m , 


<wa 


where  u 


V*' 


Since 


A-1 
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=  T  [n/^(v)  f  j  I,  (v)]  ^  (h) 

where  N2(v)  and  J2(v)  are  the  Neumann  and  Bessel  functions  of 
the  second  order,  respectively,  (A5)  may  be  written 

e  dbJ  (ai) 


-  «« 


v')  C04  Z.V  -t  2  V  f/M  iV^ 


when  t>o.  Performing  some  of  the  Integration,  (A7)  Is  written 


where 


(*  • 


and 


^<0  rtV*" 

^ (?) '/ v^[ViViv (v)^'  e 

Finally,  since  aJ  (o,t)  -  o  for  t<o. 


X  C")*)  ~  J" Ax  . 


(/If) 


(tio) 


A-2 
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APPENDIX  B 


1.  By  substituting 


H^(o,u))  =  f  dk  (ii) 

In  (25),  the  Inverse  Fourier  transform  may,  upon  changing  the 
order  of  Integration,  be  written 


(6,f)  =y  8^  aA  , 


(gj.) 


where 

B? 


(Jid><»d*)^,  and  K]^(yd)  and  K2(l^d)  are  the  modified  Bessel 
functions  of  the  second  kind  of  the  first  and  second  order, 
respectively.  Noting  that  arg  t  Ju)  ■  ±  (B3)  may  be 

written  ^ 


When  t<o,  the  impulse  response,  B™  (o,t).  Is  equal  to  zero 
since  it  Is  assumed  that  Hy  (o,t)  •>  o  for  t  A  o.  Integrating 
(B2)  by  parts.  It  follows  that 


B-l 
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where 


e  tfu. 


(B<) 


When  t  >o,  u)  In  (B6)  is  repleced  by  the  complex  variable 
B  ■  M>  -f  J  u  ,  end  the  integretlbn  Is  performed  around  the  con¬ 
tour  shown  In  Figure  9.  When  e o  and  , 

-le (6t) 


where  u  ■  — 2^ 


/<.<r 


W  * 


Since 


and 


KiOv).  -  f  [j  (V  )-i  V,  (v)  J 

?^W.(v)+j  J,(v)j  ^ 


(u) 


(if) 


where  J2(V)  and  N2(V)  are  the  Bessel  and  Neumann  functions  of 
the  second  order,  respectively,  (B7)  may  be  written 


V,(v)  -  Co«  V  Jl(v)j-v'’^inv  V,  (v)  -  as  V 


B-2 
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^ ‘  ^ - 1. (<* ) 


where 


(Tr|3)H.  ifi 

ft  \ 


,  erf  (x) 


and 


* 

v-  pin  V  w,  (v)  -  C»s V  I  (d)]  e'  <<v . 


(0.l) 


^'(**^  *  W^(v)-fasv  J^(v)je  ^  d-0  1  (gi>-) 


(Su) 
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